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To characterize the novel quantum phase transhion for a hybrid system consisting of an array 
of coupled cavities and two-level atoms doped in each cavity, we study the atomic entanglement 
and photonic visibility in comparison with the quantum fluctuation of total excitations. Analytical 
and numerical simulation results show the happen of quantum critical phenomenon similar to the 
Mott insulator to superfluid transition. Here, the contour lines respectively representing the atomic 
entanglement, photonic visibility and excitation variance in the phase diagram are consistent in the 
vicinity of the non-analytic locus of atomic concurrences. 
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Introduction: It is crucial in the modern theory of sec- 
ond order phase transitions to introduce 'order param- 
eter', whose non-vanishing average value characterizes 
one or more phases and usually breaks a symmetry of 
the Hamiltonian. But for 'quantum' phase transitions 
on the behavior of matter near zero temperature it 
is more subtle to use the traditional Ginzburg-Landau- 
Wilson paradigm, since in some cases the natural descrip- 
tion of the quantum criticality is not based on the order 
parameter [3]. 

Actually, for some systems with complicated structures 
it might be difficult to chose an appropriate order param- 
eter to correctly characterize the emergent phenomena. 
The purpose of this paper with a specific example is to 
demonstrate that, though we can not make sure what is 
the appropriate order parameter for a photon-atom hy- 
brid system, some physical observable quantities can be 
used to witness its quantum critical phenomenon. 

The hybrid system we consider is a coupled waveguide 
resonator array (CWRA) where each cavity is doped with 
a two-level atom (see Fig. 1). This hybrid architecture 
was suggested as a quantum coherent device to transfer 
and store quantum information as well as to create the 
laser-like output As for the quantum phase tran- 

sition, it is observed that such a doped CWRA can simu- 
late the Mott like transition of light from "the Mott insu- 
lator (MI) to superfluid (SF)" [6| since a doped atom can 
induce the effective photon-photon interaction in each 
cavity. Together with the inter-cavity hopping of local- 
ized phonons, this nonlinear photon-photon coupling can 
result in the Bose-Hubbard model for Mott phase tran- 
sition Recent experiments Q using cold atoms in an 
optical lattice have clearly demonstrated the quantum 
phase transition predicted by the Bose-Hubbard model. 
Actually the Bose-Hubbard theory of Mott phase transi- 
tion for cold atoms [oj is also based on the assumption of 
the order parameter, the average of the annihilation op- 
erator of boson in each site. In mean field approach, the 
average of the annihilation operator of boson is usually 
employed, while "number variance" is used in many other 




FIG. 1: Schematic setup of a cavity array with each one con- 
taining a two-level atom. Photons of mode uja can tunnel 
between adjacent cavities with hopping integral t and couple 
to the atoms with strength g. This atom-photon lattice is ex- 
pected to simulate Mott insulator and superfluid transition. 



methods to discriminate Mott insulating and superfiuid 
phases [13, 

For the hybrid system, the study of quantum phase 
transition of light in Ref. [6| still assumes the same order 
parameter in terms of photons, while more strict assump- 
tion of order parameter 1^, 11 1 was implicitly made in 
terms of the number of the polariton, which is a mixture 
of photon and atom. Here, we will not adopt such local 
order parameters of quasi-particle as direct characteri- 
zations of quantum phase transition, but pay attention 
to some observable quantities to characterize the critical 
phenomenon of a hybrid system. To this end we make 
use of atomic entanglement from the point view of quan- 
tum information, as well as photon visibility in terms of 
quantum optical theory. We remark that, as a quantum 
nonlocal property, the quantum entanglement plays an 
important role in the study of quantum phase transitions 



[iJlillS, 14, 15, 3]. We examine the signature of the Mott 
insulator to superfluid transition, the excitation number 
variance, and other two observable quantities, concur- 
rence between two atoms and the visibility of photons, 
in lattice atom-photon hybrid systems of small size by an- 
alytical and numerical methods respectively. Our results 
reveal nontrivial connections among the three quantities 
in such an intriguing way: contour lines of three quan- 
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titles in the phase diagram are approximately consistent 
with each other when the non-analyticity of concurrences 
occurs. It firmly shows that such three quantities are sig- 
natures for the MI to SF transition in such a atom-photon 
hybrid system. 

Model setup and the quasi- excitation fluctuation. We 
consider an array of N coupled cavities with each one 
containing a single two-level atom [ESS, [HI. The pho- 
ton transmission of cavity-to-cavity occurs as a hopping 
mechanism if there were the interaction between atom 
and cavity mode. Such hybrid system can be imple- 
mented with the defect array in photonic crystal [17[ 
or Josephson junction array in cavity [3|. The Hamilto- 
nian of a hybrid system, or a lattice atom-photon system, 
H = Hf ree + Hint + Hhop is decomposed as three parts, 
free Hamiltonians of light and atom, 

N N 

Hfree=uJa'^a\a,+uJb'^\e)^{e\, (1) 

1=1 1=1 

the cavity-modc-atom interaction in the ith defect 

N 

J?mt=5E(«llff)»(eKH.c.), (2) 

with strength g and the photon hopping between NN 
defects 

N 

Hhop = + H-c.) , (3) 

2 = 1 

with hopping integral constant t for the tunneling be- 
tween adjacent cavities. Here, \g)^ (|e)-) denotes the 
ground (excited) state of the atom placed at ith cav- 
ity; a\ and Oi are the creation and annihilation operators 
of a photon at defect i. Obviously the total excitation 
number 

■^=E-^' = E + + (4) 

i=l 2=1 ^ 

is conserved quantity for the Hamiltonian H, i.e., 
[H,N] = 0, where 25fje), = |e), and 25f \g), = - \g)^. 

It can be seen that M is just the single excitation num- 
ber of the polaritons. It is well known that the con- 
ventional Mott insulator to superfluid phase transition 
occurs in a Bose-Hubbard model. Here, when the repul- 
sive interaction between bosons is large enough in the 
Mott phase, the number fluctuation would become ener- 
getically unfavorable, forcing the system into a number 
state and exhibiting vanishing particle number fluctua- 
tion. In the superfluid regime, atoms are delocalized with 
non-vanishing particle number fluctuation. As for the 
present hybrid system, the fundamental excitations are 
polaritons [lH and the mechanism of the Mott transition 
is due to the effect of photon blockade. Since the pho- 
ton number is not conserved in such system, the photon 



number fluctuation An^ = A(a|^ai) is not appropriate to 
characterize the superfluid phase as that for a pure Bose- 
Hubbard model. This is because An^ does not vanish 
even in the Mott insulator regime due to the couplings 
between photons and atoms. Hereafter, we define the 
variance Ayl by (Av4)2 = {{Af) - {Af. Therefore, one 
can take the excitation number fluctuation per site AA/i 
as an order parameter to characterize the Mott transi- 
tion. In the large detuning limit (5 = lJq — cjb ^ 0, all 
atoms are in excited states, which is perfectly number 
squeezed states, i.e., AA/i = for all sites. In the other 
limit i5 <C 0, all atoms are in ground states. Obviously 
two-atom concurrence vanishes and the density fluctua- 
tion becomes AAj = \J {a\a\a.,ai) = ^/{N - 1)/A^ ~ 1 
since N = N in this case. 

Atomic entanglement characterized by concurrence. 
Intuitively, two atoms in two adjacent cavities should en- 
tangle with each other due to the hopping of phonon from 
one cavity to another. Now we try to describe this kind of 
atomic entanglement induced by coupled photons. Obvi- 
ously, if the photon is in quantum phase transition, the 
critical photon induced atomic entanglers can character- 
ize this critical behavior. 

We express the concurrence characterizing quantum 
entanglement in terms of observable quantities such as 
correlation functions. The complete basis vectors of the 
total system are denoted by 

N 

\{nj,Sj}) = \ni, ..,nN;si, ..,sn) = J| \nj) (g) \sj) (5) 

i=i 

where \nj) is the Fock state of photon and \sj) ~ \g)^, |e),j 
for Sj = 0, 1 respectively. The fact that J\f is conserved 
can be reflected by the matrix element vanishing of the 
density operator p — p{H) on the above basis for any 
state of the hybrid system, that is. 

The functional p{H) of the Hamiltonian may be a ground 
state or thermal equilibrium states. The reduced density 
matrix p^^^^ = TvpTr^ Ar[p(^)] for two atomic quasi- 
spins, e.g., si and S2 are obtained as 

r„(i2)i , , _ n^'-"'^"'^"^--'" S\S^(h- ~ <i')] 

If Js'jS^jSiSa — Pnj,sxS2S3...SN'^[/ ,\^J 

[llj;S3...SN\ 

= d{Sl+ S2- S-^- S2) Pn],s\As3... 

[nj;S3...sjv] 

(7) 

by tracing over all photon variables (with Trp) and 
atomic variables except for si and S2- The corresponding 
reduced density matrix for two atoms i and j is of the 
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form 



^(y) 



f u+ Q \ 



V / 



z*- 
Zij wjj 



(8) 



the concurrence Ca = 



According to Refs. 
max{0, Ai — A2 — A3 — A4} shared between two atoms i 
and j is obtained in terms of the the square roots {A^} 
(Ai = max{Ai}) of eigenvalues of the non-Hermitian ma- 
trix pp. Here p = {ay ® Gy) p* {uy ® Gy). Using the 
observable quantities, the quantum correlation Zij = 
(^1 SfS- IV'), u% = (1/2 ± Sf) (1/2 ± SI) 1^), the 
concurrence is rewritten as a computable form 



Cij = 2 max(0, 



(9) 



We note that this formula for the concurrence of two 
quasi-spin in a hybrid system is the same as that for 
pure spin-1/2 system [3 The non-analyticity of 

concurrence arises from the abrupt switch of the sign of 

quantity \zij\ — \J~u^tjU^j and can be used to determine 
quantum phase transitions. 

Photon visibility in hybrid system. Similar to the tran- 
sition of superfluid to Mott insulator in Bose-Hubbard 
model j20| , two phases of the atom-photon hybrid system 
can also be delimited through the quantum coherence of 
the ground state. In Mott insulating phase, the quantum 
coherence of photons is completely destroyed due to the 
photon blockade. In superfluid phase limit, the quan- 
tum coherence of photons gets its maximum. Therefore 
the quantum coherence of photons can be employed to 
indicate phases, which is characterized by a observable 
quantity, the visibility of 'interference fringes' 



(10) 



Here, Knax and Vmin are the maximum and minimum of 
the photon number distribution of ground state in the k 
space 



V{k) 



(11) 



In the strong photon blockade limit, V = while in the 
superfluid limit, V = 1. Comparing with the local quan- 
tity, the photon number fluctuation, the visibility is more 
appropriate to discriminate two phases. 

Characterizing quantum criticality. For the lattice 
atom-photon system, we now consider the connections 
among three quantities AAfi , Cij , and V around the crit- 
ical point. 

We start with an extreme case <? = 0. At zero tem- 
perature, the physics of the lattice atom-photon model 
can be described in two regimes separated by the bound- 
ary 6 = 2t. In the region 6 > 2t {S < 2t), the ground 



state is in typical Mott insulating (superfluid) phase with 
AM = i^iN- l)/N), V = (1), and = (0), 
respectively. At line S — 2t, the model admits multi- 
fold degenerate ground states with energy e^^^ = Nujb, 
and the excitation number fluctuation and visibility ex- 
perience a big jump, while the concurrence between two 
atoms is 'uncertainty' due to the energy-level crossing. 

When the atom-photon interaction g is switched on, it 
becomes avoided level crossing. This fact will result in 
the quantum fluctuation driving the transition from Mott 
insulator to superfluid phase, which corresponds to the 
non-vanishing concurrence between atoms. To illustrate 
this mechanism quantitatively, we just switch on atom- 
photon couplings in cavities i and j and leave all other 
coupling to be zero. For very small g, the un-perturbable 
ground states 



61) = |^).=o n 1.9)/ ' l'^2) = \N- l),.o IQ), 



103) = |iV-l>fe=o|e),|G,), 

104) = |iV-2)fe=o|e)Je),|G,,,), 



(12) 



are degenerate, where \n)k denotes the photon Fock 
state in k space and \Gij^,,) — Yii^i j If); denotes the 
atomic state of all atoms except I — i,j, .... Up to the 
first order perturbation with energy correction e^^"^ — 
-g^/2N(TTW) = -gVNrr^ where /3 = ^/(W^^T)/N, 
the perturbed ground state is 



|^3)=rK|0i)+/3|04))- 2(l</'2) + |03)). (13) 
The corresponding concurrence can be calculated as 



" 2(l + /32)- 



(14) 



As 5 being apart from the degenerate point, the concur- 
rence Cij decreases due to the energy competition of two 
phases. Therefore, this heuristic analysis has shown the 
simple relation among concurrence, visibility, and exci- 
tation number fluctuation around quantum phase transi- 
tion critical point: the excitation number fluctuation and 
visibility both exhibit an abrupt drop while the concur- 
rence has a sharp maximum. It can be predicted that as 
g increases, changes of the three quantities will be slow 
due to the strongly coupling between atoms and photons. 
In the following, it will be investigated for small system 
in wide range of parameters by numerical simulations. 

We investigate three quantities in a small size sys- 
tem by exact diagonalization method. For open chain 
cavity array system, the visibility V can be calculated 



where 



by S{k) = 2/{N + 1) sin(fcz) sm{kj) 

k — mi/ {N + 1), n e [1, A^], while the concurrence and 
excitation number fluctuation can be characterized as av- 
erage concurrence C = {^/N)J2i<:j ^'^'^ average ex- 
citation number fluctuation 'KM — l/A^X^i KN'i. In Fig. 
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FIG. 2: (Color online) Contours of three quantities, AAf 
(dark lines in (a,c)), V (dark lines in (b,d)), and C (color 
maps in (a-d) obtained by exact diagonalization for 2- (a,b), 
4- (c,d) cavity systems. Red lines in (a-d) denote closer con- 
tour lines of V and A7V to the non- analytical curve of C . 
It shows that contour lines of three quantities are consistent 
in the vicinity of the non-analytic locus, at which the non- 
analyticity of concurrence occurs. 



2, contours of three quantities obtained by exact diag- 
onalization are plotted in the 5/g-t/g plane for 2- (Fig. 
2(a, b)), 4- (Fig. 2(c, d)) cavity systems. Contours of ex- 
citation number fluctuation AA^ (dark lines in Fig. 2(a, 

c) ) and visibility of photons V (dark lines in Fig. 2(b, 

d) ) are compared with the concurrence C (color maps in 
Fig. 2 (a-d)) as functions of the scaled detuning 5/ g and 
photon hopping integral t/g. We see that contour lines of 
three quantities are consistent in the vicinity of the locus 
at which the non-analyticity of concurrence occurs. The 
non-analytic locus in 5/g-t/g plane is defined by the equa- 
tion \zij \ — ^Jufj'^ij — 0- Red lines in Fig. 2(a-d) denote 

closer contour lines of V and AA/" to the non-analytical 
curve of C. It also shows that the visibility and excita- 
tion number fluctuation start to drop at the non-analytic 
locus of concurrence. There is a slight difference between 
profiles of 2 and 4-cavity systems. The red contour line 
of V in 4-cavity system is closer to the non-analytic lo- 
cus of C than that in 2-cavity system. It indicates that 
contour lines of three quantities will cover at the vicinity 
the non- analytic locus of C in thermodynamics limit. 

Summary: In summary, we have investigated the Mott 
insular to superfluid transition in a hybrid system con- 
sisting of an array of coupled cavities doped with two 
level atoms. We investigate two non-local observable 
quantities, the concurrence between atoms and the vis- 
ibility of photons, in comparison with the local order 



parameter, excitation number fluctuation, for the Mott 
insulator to superfluid transition. It can be observed 
form analytical and numerical simulation results that 
the atomic entanglement and photonic visibility in the 
phase diagram indeed can reflect the quantum critical 
phenomenon signatured by the total excitation variance. 
In principle, such non-local observable quantities of the 
hybrid system can be used to detect the critical point in 
experiment. 
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